Quantum Szilard Engine 



O 

(N 

C 

00 



43 

i ' 

-i— > 

£ 



> 

O 
O 



Sang Wook Kim 1 ' 2 , Takahiro Sagawa 2 , Simonc De Liberate 2 , and Masahito Ucda 2 
1 Department of Physics Education, Pusan National University, Busan 609-735, Korea and 
2 Department of Physics, University of Tokyo, Tokyo 113-0033, Japan 
(Dated: January 11, 2011) 

The Szilard engine (SZE) is the quintessence of Maxwell's demon, which can extract the work from 
a heat bath by utilizing information. We present the first complete quantum analysis of the SZE, 
and derive an analytic expression of the quantum-mechanical work performed by a quantum SZE 
containing an arbitrary number of molecules, where it is crucial to regard the process of insertion 
or removal of a wall as a legitimate thermodynamic process. We find that more (less) work can be 
extracted from the bosonic (fermionic) SZE due to the indistinguishability of identical particles. 

PACS numbers: 03.67.-a,05.30.-d,89.70.Cf,05.70.-a 



Maxwell's demon is a hypothetical being of intelligence 
that was conceived to illuminate possible limitations of 
the second law of thermodynamics [H, 0|. Leo Szilard 
conducted a classical analysis of the demon, considering 
an idealized heat engine with a one-molecule gas, and 
directly associated the information acquired by measure- 
ment with a physical entropy to save the second law Q. 
The basic working principle of the Szilard engine (SZE) 
is schematically illustrated in Fig. [T] If one acquires the 
information concerning which side the molecule is in after 
dividing the box, the information can be utilized to ex- 
tract work, e.g., via an isothermal expansion. The crucial 
question here is how this cyclic thermodynamic process 
is compatible with the second law; the entropy of ks In 2 
(fcs is the Boltzmann constant) that the engine acquires 
during the isothermal process is not returned to the reser- 
voir but seems to be accumulated inside the engine. Now 
it is widely accepted that the measurement process in- 
cluding erasure or reset of demon's memory requires the 
minimum energy cost of at least fee In 2, associated with 
the entropy decrease of the engine, and that it saves the 
second law 



Although the SZE deals with a microscopic object, 
namely an engine with a single molecule, its fully quan- 
tum analysis has not yet been conducted except for the 
measurement process HQ. In this Letter we present 
the first complete quantum analysis of the SZE. In the 
previous literature, it takes for granted that insertion or 
removal of the wall costs no energy. This assumption is 
justified in classical mechanics but not in quantum me- 
chanics because the insertion or removal of the wall 
alters the boundary condition that affects the eigenspcc- 
trum of the system. As shown below, a careful analysis of 
this process leads to a concise analytic expression of the 
total net work performed by the quantum SZE. If more 
than one particle is present in the SZE, we encounter the 
issue of indistinguishability of quantum identical parti- 
cles. Indeed, how much work is extracted from the quan- 
tum SZE strongly depends crucially on whether it con- 
sists of either bosons or fermions. We also show that the 
crossover from indistinguishability to distinguishability 
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FIG. 1: (Color online) Schematic diagram of the thermo- 
dynamic processes of the classical SZE. Initially a single 
molecule is prepared in an isolated box. (A) A wall depicted 
as a vertical gray bar is inserted to split the box into to two 
parts. The molecule is represented by the dotted circles to 
indicate that at this stage we do not know in which box the 
molecule is. (B) By the measurement, we find where the 
molecule is. (C) A load is attached to the wall to extract a 
work via an isothermal expansion at a constant temperature 
T. 



gradually occurs as the temperature increases. We as- 
sume that the measurement is performed perfectly. The 
case of imperfect measurement is discussed in terms of 
mutual information in Ref. Q. 

To define the thermodynamic work in quantum me- 
chanics, let us consider a closed system described as 
Hip n = E n ip n , where H, ijj n and E n arc the Hamiltonian 
of the system, its nth eigenstate and eigenenergy, respec- 
tively. The internal energy U of the system is given as 
U = E n P n , where P n is the mean occupation number 
of the nth eigenstate. In equilibrium P n obeys the canon- 
ical distribution. From the derivative of U, one obtains 
dU = ^2 n {E n dP n + P n dE n ). Analogous to the classi- 
cal thermodynamic first law, TdS = dU + dW, where S 
and W are the entropy and work done by the system, 
respectively, the quantum thermodynamic work (QTW) 



2 



can be identified as dW = -J2 n P n dE n EH El- Note 
that J2 n E n dP n should be associated with TdS since the 
entropy S is defined as S = — k B YJ n Pn m Pn- 

Although the process of inserting a wall is accompanied 
by neither heat nor work in the classical SZE, it is not the 
case with the quantum SZE. This process can be modeled 
as that of increasing the height of the potential barrier. In 
quantum mechanics, energy levels then vary, contributing 
to the QTW. This process can be performed isothermally 
so that the temperature is kept constant during the whole 
process in conformity with the original spirit of the SZE. 
If the insertion is performed in an adiabatic process de- 
fined as dQ = ^2 n E n dP n = 0, one can easily show that 
the temperature is either changed or not well-defined at 
the end. The former is obvious considering the classical 
thermodynamic adiabatic process. If the wall is inserted 
in a quantum adiabatic manner, dP n = is always sat- 
isfied. Given that the temperature is defined from the 
ratio of probabilities as P n /P m = e ~ ( - En ~ Em ^ kB ' F \ it is 
well-defined only if all energy differences are changed by 
the same ratio [T3|. However, this cannot be achieved 
in the SZE because each energy level shifts in a different 
manner 
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To describe the quantum SZE it is indeed sufficient 
to know only the isothermal process for a whole cycle. 
If the external parameter X is varied from Xi to X 2 
isothermally, the QTW is obtained as 



W = k B rY^ 



x 2 



dlnZdEr, 



dX 



Xl dE n dX 
k B T[\nZ(X 2 )-\nZ(X 1 )}, 



(1) 
(2) 



-0E n 



is the partition function with ft 



where Z 
l/k B T. 

It is worth mentioning that the isothermal process in- 
duces thcrmalization of the molecule with the reservoir at 
every moment, which destroys all the coherence among 
energy levels. Therefore, it is not necessary to describe 
the dynamics of our system in terms of the full density 
matrix; it is sufficient to know its diagonal part, i.e. the 
probabilities P n . However, this thcrmalization has noth- 
ing to do with the measurement of the location of the 
molecule since it proceeds regardless of which box the 
molecule is in. 

Let us now consider a quantum SZE in a general situ- 
ation. As shown in Fig. [2] the whole thermodynamic cy- 
cle consists of four processes, namely insertion, measure- 
ment, expansion and removal among four distinct states 
(I)- (IV). N ideal identical molecules are prepared in a po- 
tential well of size L as shown in Fig. [2JI) . A wall is then 
isothermally inserted at a certain position I. The parti- 
tion function, at the moment when the wall insertion is 
completed but the measurement is not performed yet, is 
given as Z(l) = Y^ m =o ^™(0> where Z m (l) denotes the 
partition function for the case in which m particles are 
on the left side and JV-mon the right. The amount of 
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FIG. 2: Schematic diagram of the quantum SZE containing 
three molecules. (I) Three molecules are prepared in a closed 
box with size L. (II) A wall, depicted by a vertical gray bar, 
is isothermally inserted at location I. The process (I) — s> (II) 
is called 'insertion'. (Ill) The information on the number of 
molecules, m, on the left is acquired by the measurement. The 
process (II) —¥ (III) is called 'measurement'. (IV) The wall 
moves and undergoes an isothermal expansion until it reaches 
its equilibrium location denoted by ?£q. The process (III) — ► 
(IV) is called 'expansion'. Finally the wall is isothermally 
removed to complete the cycle. The process (IV) — > (I) is 
called 'removal'. 



work required for the insertion process is thus expressed 
as 



W ins = k B T[\nZ{l) -lnZ(L)]. 



(3) 



Then, the measurement is performed without any ex- 
penditure of work. The amount of work extracted 
via the subsequent isothermal expansion is given as 
W cxp = k B TJ2^ =Q f m [In Z m (l™)-]nZ m {l)] , where 
f m = Z m (l)/Z(l) represents the probability of having 
m particles on the left at the measurement. The wall 
moves until it reaches an equilibrium position deter- 
mined by the force balance, F left + F rlght = 0, where the 
generalized force F is defined as P n (dE n /dX), as il- 
lustrated in Fig. HJIV). We note that is not simply 
(m/N)L unlike classical ideal gases. 

The wall is then finally removed isothermally. In real- 
ity the wall is not impenetrable, and has a finite poten- 
tial height, namely X^. During the expansion process, 
Xoo is assumed to be large enough to satisfy r t ^> t , 
where r t and r are a tunneling time between the two 
sides and an operational time of thermodynamic pro- 
cesses, respectively, ensuring that m is well defined. Dur- 
ing the wall removal, however, r t gradually decreases and 
becomes comparable with r for a certain strength, Xq, 
where any eigenstate is delocalized over both sides due 
to tunneling. It implies that the partition function is 
given by Z(J£) = £^ = o^(Q) rather than Z m (l%). 
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TABLE I: Total work measured in units of ksT of the quan- 
tum SZE with I = L/2 containing two bosons or two fermions 
at the low and the high temperature limits (see pH ] for de- 
tailed derivation). 
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The integral ([T]) for each to can be split into two parts; 
f*° [d\nZ m (l™^/dX]dX and J° o [dlnZ(l™)/dX]dX. It 
is then shown that the former vanishes as far as the quasi- 
static process, r — 5- oo (i.e. A , X^ — > oo), is concerned. 
This leads us to 



k B Tj2f m [lnZ(L)-lnZ(l™)] 



(4) 



m=0 



Note that the summation over to must be made in Z(l) of 
Eq. ([3]) for the insertion process irrespective of tunneling 
since no measurement has been performed yet. 

Combining all the contributions the total work per- 
formed by the engine during a single cycle is given by 



N 

-k^T ^2 f m In 

rn— 



(5) 

where f* m = Z m (l™)/Z(l™). Equation © has a clear 
information-theoretic interpretation in the context of 
the so-called relative entropy or Kullback-Leibler 
divergence even though f m is not normalized, namely 
J2 m fm 1 ■ K nas been shown that the average diss- 
pative work upon bringing a system from one equilib- 
rium state at a temperature T into another one at the 
same temperature is given by the relative entropy of the 
phase space distributions between forward and backward 
processes With filtering or feedback control like mea- 
surement processes of the SZE that determine to, the 



work is represented as a form equivalent to Eq. ([5]) 17 [. 

It is straightforward to apply Eq. (0 to the original 
SZE consisting of a single molecule of mass M. For 
simplicity let us consider an infinite potential well of 
size L, and I = L/2. One finds /q = /jf = 1 since 
in these cases the wall reaches the end of the box so 
that Z(r c n q ) = Z m (l™) (to = 0,1) is satisfied. Note that 
f m = 1 is always true for to = and N. Together with 
fo = fi = 1/2, wc obtain W to t = fc B Tln2, implying 
the work performed by the quantum SZE is equivalent 
to that of the classical SZE. However, consideration of 
individual processes reveals an important distinction be- 
tween the classical and quantum SZE's. For the quantum 
SZE one obtains W ins = -A + k B T \n2, W exp = A and 
Wrem = for each process, where A = ]n[z(L)/z(L/2)], 
z(l) = En=i e ~' 3E " (l) ' and E nQ) = h 2 n 2 /(SMl 2 ) with 
h being the Planck constant. In the low-temperature 



limit, A is simply given as E\(L/2) — E\{L). If the 
insertion process were ignored in the classical SZE, the 
second law would be violated because A 3> k B T in the 
low-temperature limit. In fact, A of the expansion pro- 
cess is compensated by the work required for inserting 
the wall. In the end, a tiny difference of work between 
these two processes results in the precise classical value, 
/cB? 1 ln2. As the temperature increases, the classical re- 
sults of individual processes are recovered, i.e. W- ms — > 0, 
Wexp — > k B T In 2 and Wrem = since A approaches 
fceTln2 in this limit. 

For the quantum SZE with more than one particle dra- 
matic quantum effects come into play. Let us consider a 
two-particle SZE confined in a symmetric potential well 
and a wall at I = L/2. One also finds /q = /| = 1 
because of the same reason mentioned above. Since for 
m = 1 the wall does not move in the expansion process, 
implying I = Z^ q , one obtains /i = /*. We thus end up 
with 



-2fc B r/oln/ , 



(6) 



where fo = is used. 

To get some physical insights let us consider two lim- 
iting cases of Eq. (|6]), which is summarized in Table U 
(see [l8| for a detailed derivation). For simplicity here 
the spin of a particle is ignored. In the low-temperature 
limit only the ground state is predominantly occupied, 
so that there exists effectively only one available state 
for each side. It is clear, as shown in Fig. 3(a), that for 
bosons fo should become 1/3, i.e. W tot = &bT(2/3) hi2, 
since we consider two mutually indistinguishable bosons 
over two places. On the other hand, two fermions are pro- 
hibited to be in the same side or to occupy the same state 
due to the Pauli exclusion principle, which explains why 
the work vanishes in the low-temperature limit. How- 
ever, the higher the temperature, the larger the num- 
ber of available states. Thermal fluctuations wash out 
indistinguishability since two identical particles start to 
be distinguished by occupying different states. This is 
why in the high-temperature limit we have /o = 1/4, i.e. 
Wtot = fcBTln2, for both bosons and fermions, which 
results from allocating two distinguishable particles over 
two places as shown in Fig. 3(b). It is also shown in 
Fig. 3 that fo continuously varies from 1/3 (0) to 1/4, 
exhibiting the crossover from indistinguishability to dis- 
tinguishability for bosons (fermions) as the temperature 
increases. 

The inset of Fig. [3] clearly shows that one can extract 
more work from the bosonic SZE but less work from the 
fermionic SZE over the entire range of temperature. (See 
[lH for detailed discussions of the W to t{T).) While de- 
tails of Wtot(T) depend on the confinement potential, its 
low-temperature limits given in Table U are universal and 
have a deep physical meaning associated with the infor- 
mation content of quantum indistinguishable particles as 
mentioned above. 
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FIG. 3: (Color online) /n as a function of T for bosons (solid 
curve), fermions (dashed curve) and classical particles (dash- 
dotted line) in the case of the infinite potential well. The 
temperature is given in units of Ei(L) /k&. (a) Three pos- 
sible ways in which two identical bosons are assigned over 
two states, (b) Four possible ways in which two distinguish- 
able particles are allocated over two places. The inset shows 
Wtot /W c as a function of T. 



Finally, we briefly mention possible experimental real- 
izations of the quantum SZE. Although several experi- 



ments [191 - 4221 ] or proposals [23, [24| associated with the 



classical SZE have been presented so far, its fully quan- 
tum version has been elusive. There exist three impor- 
tant ingredients for the experimental realization: (i) con- 
trollability of the confinement potential, (ii) availabil- 
ity of a thermal heat bath to perform isothermal pro- 
cesses, and (iii) measurability of the work performed. For 
bosons, the system of trapped cold atoms may be a good 
candidate since the confinement potential can be easily 
controlled. Although such a system usually lacks a ther- 
mal heat bath, it can be immersed in a different species 
of atom trapped in a wider confinement potential so that 
they can play a role of a heat reservoir. For fermions, 
two-dimensional electron gases confined in quantum dots 
made of semiconductor heterostructures might be a can- 
didate due to its controllability of the confinement poten- 
tial and the existence of a heat reservoir of the Fermi sea 
of electrons. In principle, the work is determined once 
both E n and P n are known during the entire isothermal 
processes for both bosons and fermions. 

In summary, we have studied the quantum nature of 
the SZE. The total work performed by the quantum SZE 
is expressed as a simple analytic formula which is directly 
associated with the relative entropy in the classical limit. 
To correctly describe the quantum SZE the processes of 
inserting or removing the wall should be regarded as a 
relevant thermodynamic procedure. The quantum SZE 
consisting of more than one particle clearly shows the 
quantum nature of indistinguishable identical particles. 
We believe our finings shed light on the subtle role of 



information in quantum physics. 
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We present the derivation of Table I and the detailed discussion of Wtot{T) in the high temperature 
limit. 

PACS numbers: 



I. DERIVATION OF TABLE I 

Even though Wtot obtained from Eq. (6) depends on what the confinement potential is, in the limits T — > and 
\ T — > oo, it has a universal value under quite general conditions discussed below. We first consider the two-boson case 
for an arbitrary confinement potential V(x) with reflection symmetry, V(x) = V(—x). The wall is inserted at x = 0, 
i.e. I = L/2 separating V into two single- wells. The probability that two bosons are in the right well is given as 



o 



oo; / = ^L, (si) 



^ c"| where 

- ,(2/3) (S2) 

a ■ 

^ ' with z(/3) is the partition function of a single particle in the single- well at temperature T = (&b/?) ■ 
O" 1 To obtain universal values of /o in the high- and low-temperature limits, the key observation is that 



(N 
> 



d -4 oo (T oo) (S3) 



and 



^ i d -> 1 (T -> 0) (S4) 

1 hold only with the following two assumptions: 
' 1. The number of energy levels is infinite. 

O ! 

y—i . 2. The single-particle ground state is not degenerate. 

We first prove Eq. (|S3|) . Since z(/3) > z(2/3) holds, we have d > z(/3). On the other hand, we can show that 
, z(f3) — > oo holds with T — > oo from assumption 1. Therefore we obtain Eq. (|S3[) . Next, from assumption 2, we have 

d = {l + 2e^ l3i - El ^ Ea ^ H )/(l + e -2/3(£i--E ) _| ^ where is the nth energy level of the single-well. It immediately 

_<?3 ; leads us to Eq. ([gl]t. From Eqs. ([S3]) and ([5i]). we obtain 

/o -> 1/4 (T -> oo) (S5) 



and 



/o -> 1/3 (T ->. 0), (S6) 



which are potential-independent universal results and describe the distinguishability-indistinguishability crossover by 
varying temperature. The total work is then given as 

Wtot ->fc B Tln2 (T^oo) (S7) 

and 

Wtot (2/3)fc B Tln3 (T 0), (S8) 

which are also model-independent. 
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FIG. SI: Wtot/Wc as a function of T for bosons (solid curve) and fermions (dashed curve) obtained from Eq. (5). The dashed 
dot line represents the classical work. This is the same as the inset of Fig. 3. The inset shows Wtot — W c for bosons (solid 
curve) and fermions (dashed curve) as a function of T for an infinite potential well (thick curves) and for a harmonic potential 
(thin curves), which exhibits the deviation proportional to ±T X / 2 for the infinite potential well. (See the text for details.) For 
the harmonic potential kui — 10£a is used, where uj and Ei represent the resonant angular frequency of the harmonic oscillator 
and the energy of the ground state of the infinite potential well, respectively, and h = h/(2n). 



A similar approach can be done for the spinlcss two-fcrmion case, /o is given as 

where d is also defined as (|S2|) . After exactly the same assumptions and similar procedure done above one reaches 

/o -> 1/4 (T -> oo) (S10) 

and 

jWO(T^O). (Sll) 

Therefore, we have 

Wtot ^fc B Tln2 (T^oo) (S12) 

and 

Wtot -> (T -> 0). (S13) 



II. Wtot AS A FUNCTION OF T 



Figure [STI shows W to t/W c is greater than one for bosons and smaller than one for fermions. As the temperature 
increases, Wtot/W c approaches one for both bosons and fermions as expected. It is noted, however, that Wtot itself 



3 



does not converge to W c but exhibits rather clear deviations from W c depending on the shape of the potential as 
shown in the inset of Fig. IS1I More precisely W to t — W c diverges as iT 1 / 2 (+ for bosons, and — for fermions) for an 
infinite potential well, and keeps constant for a harmonic potential. 

These deviations can be understood by expanding Wot perturbatively. From Eqs. (|ST|) and (fS9| one obtains 

^ = z£k> (S14) 

where b = dr 1 of Eq. (|S2|) . and + for boson and — for fermions. Note that b goes to zero, i.e. fo —> 1/4, in the high 
temperature limit. The total work is then expressed as 

W tot «fc B T^ln2±^ln^ . (S15) 

It is shown that b —> T~ 1//2 for an infinite potential well implying Wtot — W c ~ ztT 1 / 2 , and b — > T^ 1 for a harmonic 
potential implying W to t — W c ~ (fc/4) ln(4/e). In fact, one finds that b approaches T~ x / a if the energy eigenvalues 
are given as E n ~ n a . For a — > 0, Wtot approaches the classical result irrespective of temperature. Since the infinite 
potential well is the steepest potential, a = 2 is the maximum possible value in practice. Therefore, at worst Wtot — W c 
exhibits iT 1 / 2 deviation. Such a deviation is regarded as quantum effect of the SZE, and may be directly observed 
in experiment. Note that b was ignored in the previous section, so that we simply had Wtot ~^ k^T lri2 as T — > oo in 
Table I. 



